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Abstract. The X-ray luminosities of the hot halo gas around simulated, Milky 
Way like disk galaxies have been determined, as a function of redshift. The X-ray 
luminosity increases significantly with redshift, in some cases as much as a factor 
30 going from z—Q to 2. Consequently, the optimal detection redshift can be JJ0.5. 
Results of fully cosmological simulations of galaxy groups and clusters, incorpo- 
rating star formation, chemical evolution with non-instantaneous recycling, metal- 
dependent radiative cooling, galactic super-winds and thermal conduction are 
presented. X-ray luminosities at z~0 are somewhat high and central entropies 
somewhat low compared to observations. This is likely a combined effect of chem- 
ical evolution and metal- dependent, radiative cooling. Central ICM abundance 
profiles are somewhat steep, and the observed level of ICM enrichment can only 
be reproduced with IMFs more top-heavy than the Salpeter IMF. In agreement 
with observations it is found that the iron in the ICM is in place already at z ~1. 
The [Si/Fe] of the ICM decreases with time, and increases slightly with radius at 
a given time. The cluster galaxies match the observed "red sequence" very well, 
and the metallicity of cluster galaxies increases with galaxy mass, as observed. 
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1. X-ray emission from the haloes of & Houck 1997; Strickland et al. 2004), NGC 

disk galaxies 2841 (Benson et al. 2000), and the Milky 

Way itself (e.g., Sidher et al. 1996; Pietz et 

Standard models of disk galaxy formation al 1998), if neglecting cases where (i) the 

require spiral galaxies to be surrounded by galactic disks show evidence of being dis- 

large reservoirs of hot gas, which should be turbed by tidal interactions, or (ii) the X- 

emitting at X-ray wavelengths, and from ray emission at off-disk distances of a few 

which gas should still be accreting on to the kpc can be attributed to processes origi- 

disk at present (e.g., White & Frenk 1991). nating in the disk such as feedback from 

The detailed properties of these haloes and star-bursts. 

their role in the formation and evolution Toft et al. (2002) extracted disk galax- 
of galactic disks remain largely unknown ies from fully cosmological hydrodynamical 
from an observational viewpoint, as the simulations and determined their present- 
haloes have so far escaped direct detection, day halo X-ray luminosities, finding con- 
Possible exceptions are NGC891 (Bregman sistency with observational upper limits 
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Fig. 1. Redshift evolution of the rest-frame 
bolomctric and 0.2-2 keV X-ray luminosi- 
ties of a simulated, Milky Way like disk 
galaxy. 



on halo emission from nearby spirals. 
Moreover, based on the predicted rates of 
mass accretion onto the disks, they sug- 
gested that spiral galaxy haloes were possi- 
bly one order of magnitude brighter in soft 
X-rays at z ~ 1 than at present (the gas 
mass cool-out rate is expected to be propor- 
tional to the halo X-ray luminosity times 
the emissivity weighted inverse tempera- 
ture of the hot halo gas — see Rasmussen 
et al. 2004 for details). Analytical mod- 
els of halo emission overpredict emission at 
z = by at least an order of magnitude 
(see Benson et al. 2000 and Toft et al. 2002) 
and thus cannot be assumed to provide a 
reliable description of the evolution in X- 
ray emission with redshift. Rasmussen et 
al. 2004 extended the work of Toft et al. by 
studying the higher-redshift properties of 
two Milky Way like disk galaxies extracted 
from the ACDM galaxy formation simula- 
tions of Sommer-Larsen, Gotz & Portinari 
2003. Figure 1 shows the redshift evolu- 
tion of the (rest-frame) bolometric and 0.2 - 
2 keV (X-ray) luminosity of the halo of a 
Milky Way like disk galaxy. The simulation 
was run with a radiative cooling function 
corresponding to 1/3 solar abundance. The 
increase in luminosity is about a factor of 
30 going from z—0 to 2. Given this, it is 
obviously of interest to determine the opti- 



Fig. 2. 0.3-2 keV surface brightness of the 
galaxy from Fig. 1, averaged over 40 kpc 
wide slabs oriented parallel the disk, as a 
function of vertical distance |z| to the disk. 
The galaxy is viewed edge-on; for each ver- 
tical bin, the values above and below the 
disk have been added. 



mal redshift for X-ray halo detection, tak- 
ing redshifting of the photons, cosmolog- 
ical (1 + z) 4 surface brightness dimming, 
and Galactic HI absorption into account. 
Rasmussen et al. find that it varies from 
galaxy to galaxy, but in some cases the op- 
timal detection redshift can be <;0.5. Figure 
2 shows for one such case (the galaxy from 
Fig. 1) the predicted, edge-on 0.3-2 keV 
surface brightness profiles at various red- 
shifts. An observer at the current epoch 
would measure the largest surface bright- 
ness of this galaxy at |z|^10 kpc, when it 
was at z ^0.8. 

2. Simulating galaxy groups and 
clusters — the code and the 
simulations 

Only fairly recently has it been possi- 
ble to carry out fully cosmological gas- 
dynamical/N-body simulations of the for- 
mation and evolution of galaxy clusters at 
a level of numerical resolution and physi- 
cal sophistication that the cool-out of gas, 
star-formation and chemical evolution in, 
and outflows from, individual cluster galax- 
ies can be modelled to, at least some, degree 
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of realism (e.g., Valdarnini 2003; Tornatore 
et al. 2004). 

We have recently engaged in under- 
taking such simulations of galaxy groups 
and clusters. Scientific objectives in- 
clude reproducing the observed luminosity- 
temperature relation for groups and clus- 
ters and its evolution with redshift; sur- 
face brightness and entropy profiles; and 
the high iron abundance of ~l/3 solar in 
the intracluster medium (ICM). The met- 
als have to be first produced in the galaxies 
and subsequently expelled and mixed into 
the ICM; since the mass of stars in galax- 
ies is only about 10-20% of the mass of 
hot ICM gas, this requires high efficiency in 
metal production and expulsion. Moreover 
one should reproduce observed ICM abun- 
dance profiles and a-to-iron ratios, "cold 
fractions" , and group and cluster galaxy lu- 
minosity functions and "red sequences" . 

The code used for the simulations 
is a significantly improved version of 
the TreeSPH code used previously for 
galaxy formation simulations (Sommer- 
Larsen, Gotz and Portinari 2003): Firstly, 
we have adopted the "conservative" en- 
tropy eq. solving scheme suggested by 
Springel & Hernquist (2002). Secondly, 
non-instantaneous chemical evolution trac- 
ing 9 elements (He, C, N, O, Mg, Si, S, Ca 
and Fe) has been incorporated in the code 
following Lia, Portinari & Carraro (2002). 
Atomic radiative cooling depending both 
on the metal abundance of the gas and the 
meta-galactic UV field is invoked together 
with simplified radiative transfer. Thirdly, 
continuing, star-burst and/or AGN driven 
galactic super-winds are incorporated in 
the simulations. This is required to get the 
abundance of the ICM to the observed level 
of about 1/3 solar in iron. The strength 
of the super-winds is modelled through a 
free parameter / W j n d which determines how 
large a fraction of the stars born partake 
in super-wind driving star formation. We 
find that in order to get an iron abundance 
in the ICM comparable to observations, 
/wind <; 0.5 and at the same time a fairly 
top-heavy Initial Mass Function (IMF) has 
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Fig. 3. Lx — Tx relation at z=0. Data 
points are marked by open symbols, simula- 
tions by filled symbols and crosses. Details 
about the simulations are given in the text. 



to be used — see section 3. Finally, thermal 
conduction was implemented in the code 
following Cleary & Monaghan (1999). 

Six groups and clusters (T=l-6 keV) 
were selected from a cosmological, L=15Q 
ft _1 Mpc, "standard" ACDM N-body sim- 
ulation, and re-simulated with baryonic 
physics included, using the code described 
above, for a suite of physical parameters. 
The wind strength parameter / w i n( j was 
set to 0.7 and Salpeter (Sal) and Arimoto- 
Yoshii (AY) IMFs were adopted. Some sim- 
ulations were run with galactic super-wind 
feedback two times (SWx2FB) and four 
times (SWx4FB) as energetic as is available 
from supernovae; the remaining amount of 
energy is assumed to come from AGN ac- 
tivity. Others were run with an AY IMF, 
strong feedback (SW) and additional pre- 
heating at z=3 (AY-SW + preh.), as dis- 
cussed by Tornatore et al. (2003). One se- 
ries of simulations was run with a Salpeter 
IMF and only early (z^4), strong feedback, 
as in the galaxy formation simulations of 
Sommer-Larsen, Gotz and Portinari (2003) 
(as this results in overall fairly weak feed- 
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Fig. 4. ICM iron abundance profiles for 
various simulations of the "Virgo" clus- 
ter, and one (AY-SW) of the "Mini-Coma" 
cluster (long-dashed line). Also shown are 
data points from De Grandi et al. (2003) for 
cool-core clusters (open circles with error- 
bars) and non cool-core clusters (filled cir- 
cles with errorbars). The appropriate com- 
parison of the models is to the cool-core 
clusters. 

back we denote this: Sal-WFB). Another 
was run with thermal conduction included 
(AY-SW + Therm. Cond.), assuming a con- 
ductivity of 1/3 of the Spitzer value (e.g., 
Jubelgas et al. 2004). Finally, for compar- 
ison with previous work, a series of simu- 
lations was run without radiative cooling, 
star formation and thermal conduction. 

3. Simulating galaxy groups and 
clusters — results 

In Figure 3 is shown the z=0 Lx — T rela- 
tion for the 6 groups and clusters (T=l- 
6 keV), re-simulated with various physi- 
cal parameters, as explained in the previ- 
ous section. Also shown are observational 
data on the z=0 Lx — T relation from vari- 
ous sources. The predicted X-ray luminosi- 
ties are mostly somewhat high compared to 
observations, with the exception of models 
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Fig. 5. ICM iron abundance profiles for the 
AY-SW simulation of the "Virgo" cluster 
at z—0, 1 and 2 (solid lines). Also shown 
are [Si/Fe] profiles at z—0, 1 and 2 (dashed 
lines). To the right 2=2,1,0 corresponds to 
going top down for the two sets of curves 
(z=2: red; z=l: green; 2=0: blue). 



run with an Arimoto-Yoshii IMF and very 
strong feedback (AY-SWx4FB), and mod- 
els run with a Salpeter IMF and only little 
feedback (Sal-WFB), hence resulting in a 
low level of metal enrichment. The prelim- 
inary conclusion from the results of these 
simulations is that the success of previous 
generations of models in reproducing the 
Lx — T relation, the surface brightness pro- 
files, the entropy "floor" at S = T/n 2 J 3 ~ 
100-200 keV cm 2 etc. (e.g., Tornatore et 
al. 2003) likely was due to the adoption 
of a primordial radiative cooling function 
(i.e. the contribution to radiative cooling 
from metals was neglected). The analysis 
of this is ongoing, however. For most of our 
models we find central entropy values of ~ 
10-20 keV cm 2 ; this may be in agreement 
with "observed" values in cool-core clusters 
(M. Voit, this conference). 

Figure 4 shows the profiles of iron abun- 
dance in the ICM for various models, to- 
gether with recent observational data by De 
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Fig. 6. Cumulative iron mass of the ICM 
at z=0 for the Sal-SW and AY-SW simu- 
lations. Actual profiles are shown by solid 
lines, profiles inferred assuming the ob- 
served iron profiles of De Grandi & Molendi 
(2001) by dashed. Top solid line corre- 
sponds to AY-SW, bottom to Sal-SW. Top 
dashed line corresponds to Sal-SW, bottom 
to AY-SW (AY: blue; Sal: red). 



Fig. 7. Top panel: B-R vs. R Colour- 
Magnitude relation for the 94 'Mini- 
Coma' (T~6 keV) cluster galaxies at z=0, 
compared to the z~0 Red Sequence of 
Gladders et al. (1998) (straight line). 
Bottom panel: Metallicity-luminosity rela- 
tion for the "Mini-Coma" galaxies (full 
symbols: mass-wcigthed, open symbols: 
luminosity- weighted) . 



Grandi et al. (2003). All models, but one, 
were run for a T~3 keV cluster, which we 
shall denote "Virgo" in the following. One 
model refers to a larger cluster with T~6 
keV, denoted "Mini-Coma" . In general, the 
iron abundance profiles are too steep in the 
inner parts of the clusters and too low in 
the outer parts. In relation to the former 
problem, AGN induced, buoyant bubbles 
may help to transport energy and metals 
from the very central part of clusters out- 
ward, flattening the inner abundance pro- 
files. In relation to the latter problem, clear 
progress is made going from the Salpeter 
IMF to the more top-heavy Arimoto-Yoshii 
IMF — see also Portinari et al. 2004. This 
indicates that an even more top-heavy IMF 
may be required to reproduce the observed 
ICM abundance, but see also below. 

Figure 5 shows, for the "Virgo" cluster 
AY-SW simulation, the profiles of [Fe/H] 



and [Si/Fe] in the ICM, at redshifts z=0,l 
and 2. In agreement with observations the 
metals are essentially in place at z—1 (e.g., 
Tozzi et al. 2003), and it is seen that sub- 
stantial enrichment has taken place already 
at z~2. This is in agreement with the find- 
ings of Sommer-Larsen et al. (2004), that 
the cluster stars formed at typical red- 
shifts z/il2.5-3. Moreover, [Si/Fe] increases 
slightly with radius at any redshift, also in 
agreement with observations (Finoguenov 
et al. 2000). Finally, [Si/Fe] generally de- 
creases with time, as is expected due to the 
substantial, but delayed iron enrichment of 
the ICM by supernovae type la. 

Figure 6 shows the cumulative iron 
mass in the ICM for the Sal-SW and AY- 
SW simulations of the "Virgo" cluster. Also 
shown is the cumulative mass of iron as- 
suming the observed ICM iron profile of De 
Grandi & Molendi (2001) (similar to the 
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updated profiles of De Grandi et al. 2003). 
It is clear that the ICM can not be suffi- 
ciently enriched for simulations adopting a 
Salpeter IMF, whereas this can marginally 
be accomplished by using an Arimoto- 
Yoshii IMF (at least out to r~0.2-0.3r vir ). 

Concerning the properties of the cluster 
galaxies, shown in Figure 7, at z=0, is the 
colour-magnitude (B-R versus R) relation 
of the 94 galaxies in the "Mini-Coma" clus- 
ter (excluding the cD). For comparison is 
also shown the mean locus of the observed 
red sequence (solid line) from Gladders et 
al. (1998). The match to observations is 
quite satisfactory. Furthermore is shown 
the mean heavy element abundance for the 
94 galaxies versus absolute B-band mag- 
nitude. The non-zero slope of the colour- 
magnitude relation for our cluster galaxies 
mainly results from the increase of the typ- 
ical stellar metallicity with galaxy luminos- 
ity, as observed (e.g., Kodama & Arimoto 
1997). 

Full details about our results are given 
in Romeo, Portinari & Sommer-Larsen 
(2004), Romeo et al. (2004) and Sommer- 
Larsen et al. (2004). 
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